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Abstract 

A new thermometer based on fragment momentum fluctuations is presented. 
This thermometer exhibited residual contamination from the collective motion 
of the fragments along the beam axis. For this reason, the transverse direction 
has been explored. Additionally, a mass dependence was observed for this ther- 
mometer. This mass dependence may be the result of the Fermi momentum of 
nucleons or the different properties of the fragments (binding energy, spin etc..) 
which might be more sensitive to different densities and temperatures of the 
exploding fragments. We expect some of these aspects to be smaller for protons 
(and/or neutrons); consequently, the proton transverse momentum fluctuations 
were used to investigate the temperature dependence of the source. 

Keywords: nuclear reactions, temperature, phase transition, caloric curve, 
projectile fragmentation 



1. Introduction 

For many years, various thermometers have been used to expand the ex- 
perimental understanding of nuclear systems. These studies have often been 
motivated by the desire to define the proposed nuclear liquid-gas phase tran- 
sition [H, 0, 13| • A broad array of caloric curves have been obtained allowing a 
better understanding of the nuclear limiting temperature and its dependence on 
source excitation energy as well as source size [4j. Additionally, nuclear ther- 
mometry has provided a bridging connection between isoscaling [3, H, and the 
symmetry coefficient of the nuclear equation of state [|| . 
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Figure 1: Left: Raw data from AE-E with Si-Csl detectors. Right: Projected, linearized 
data for Z=4 isotopes. The x-axis is Lx , which is the linearized value, (color online) 



Several thermometers have been used extensively in nuclear studies. In many 
studies, the nuclear temperature has been obtained from energy spectra through 
moving source fitting 1Q, 11, 12| . However, this thermometer is known to exhibit 
non- thermal and collective behavior [ijjj]. Alternatively, temperature may be 
obtained through isotopic thermometers The double isotope thermometer 
has provided much data [18l . Il5l Il6l . |17| , however the temperatures derived are 
complicated by model dependent secondary decay corrections 17, 19j . Ratios of 
excited states may also be used to obtain temperature for low to mid excitation 
energies 0]. 



2. Experimental 

In this work, we seek to measure the temperature of reconstructed quasi- 
projectiles from the reactions of 86 > 78 Kr+ 64 ' 58 Ni at 35MeV/nucleon collected 
with the NIMROD-ISiS array [13,|2lj] at Texas A&M University. The NLMROD- 



ISiS array is a 4tt charged particle array housed inside the TAMU Neutron Ball. 
This experimental hardware configuration provides free neutron data in con- 
junction with isotopically resolved charged particle event reconstruction. The 
data presented in this paper was extracted from particles detected in the silicon 
detectors in rings 2-9 and the Csl detectors in rings 2-11. 

There are three possible sources of fragment charge and mass identification 
in NIMROD-ISiS. Pulse shape discrimination of the Csl fast versus slow light 
output was used for Z=l,2. The AE-E method is used for Z > 3 on Si-Si and 
Si-Csl. The particle identification was done through linearization of the raw 



data [21| to remove the non-trivial curvature due to energy deposition charac- 
teristics of the detection methods (see left panel of Fig. [1]). The linearization 
utilized lines carefully chosen to follow the strongest isotope of each element. 
The data were then straightened using a calculation of the distance between the 
data point and two closest chosen lines. 

After linearization, data were projected onto an axis, L(x), producing quasi- 
Gaussian peaks. The isotopic peaks in the projected distributions were fitted 
with Gaussian functions. The fragment charge was assigned by limits placed 
on the linearization x-axis value (Lx )■ The mass of each particle was assigned 
by determining the probability of the particle belonging to a given isotope. 
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This probability (Pa) wa s calculated by comparing the value of the the isotopic 
Gaussian functions at the Lx value of the particle 



_ G A (L X ) 

where Ga is the fitted Gaussian to the selected isotope which is compared 
to the summation over all Gaussian functions (Gj) defined for the element. 
For this data, a mass was defined only if the Pa was > 0.75. This method 
of fitting the linearized data with Gaussian functions provided the ability to 
estimate the average contamination between neighboring isotopes. The average 
contamination in the yield of a given isotope as defined here was calculated 
to be « 5% across all reaction systems and all detectors. Example linearized 
projections of Be fragments are plotted in the right panel of Fig. [TJ 

The quasi-projectile source was selected by means of several event-by-event 
cuts on the experimental data. The first cut required that the sum of the col- 
lected charged fragments (^ Z) for the event equal a minimum of Z=30. The 
fragments in an accepted event were then cut on the longitudinal velocity rel- 
ative to the largest fragment (22[. This cut varied with the fragment Z. The 
fragments retained for Z=l, Z—2, and Z>3 had to have longitudinal velocities 
within the range of ±65%, 60%, 40%, respectively, of the largest fragment longi- 
tudinal velocity. The sum of the charges of the collected and accepted fragments 
was again constrained to be in the range Z=30-34. Limits were placed on the 
deformation of the source by means of a cut on the quadrupole moment of the 



momentum distribution [23 1 



This method of source definition was compared to events and fragments 
generated by the HIPSE-SIMON code [Hj]. The initial J2 z cut eliminates a 
significant number of incomplete events and events from other sources. This 
cut, however, was not sufficient to eliminate all fragments resulting from mid- 
velocity or pre-equillibrium sources. The fragment velocity cut successfully 
helped minimize fragments from non-projectile-like sources. The final cut on 
the deformation provided a reasonable level of isotropy in the selected events. 

After the quasi-projectile source was defined, the event excitation energy 
was calculated using: 

Kource = E K CP(*) + M « (*») " Q B 9' (2) 
i 

The excitation energy (E* ource ) was defined as the sum over the accepted par- 
ticles center of mass kinetic energy (K cp and K„) minus the reaction Q-value. 
The average kinetic energy of the neutrons was calculated as the proton average 



kinetic energy with a correction for the Coulomb barrier energy 25 1. 

Proper calculation of the excitation energy should include the neutron and 
gamma energy corrections. The Neutron Ball provides an experimental neutron 
multiplicity per event that must be corrected for background, neutron source, 
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and efficiency. The multiplicity of free neutrons was extracted from the experi- 
mental quasiprojectile (Q p ) data using: 



M QP 




(3) 



(tQP + (-7/.6) 



The multiplicity of neutrons assigned to the projectile source (Mqp) was calcu- 
lated from the background (Mws S ) corrected experimental neutron multiplicity 
(M exp ). and corrected using the relative efficiency of the Neutron Ball for de- 
tecting free neutrons emitted from a quasi-projectile (egp) and quasi-target 
(cqt) sources for this reaction. The free neutron correction also accounted for 
the total neutron contributions from both the target (Nt) and projectile (Np) 
nuclei. The efficiencies were extracted from tagged neutrons generated by the 
HIPSE-SIMON code and the GEANT3 [H simulation of the detector 20]. Sim- 
ulation and filtering provided £qp=0.65 and eQp=0.40. The Neutron Ball was 
calibrated with a 252 Cf source to 70% efficiency for this experiment. The cor- 
rection factor in the denominator accounts for the experimental calibration of 
the Neutron Ball to 70% efficiency rather than the 60% efficiency obtained from 
the GEANT3 simulation of the detector. Event by event, the Mqp provided an 
estimate of the free neutrons emitted from the quasi-projectile. 

The energy carried away by gammas is generally considered to be small [3, 
[27] | . In an effort to account for gamma energy, the ISiS collaboration calculated 
E 1 = lMeV(Mz>3) [13]. However, this energy has not been constrained ex- 
perimentally [281 ]. Using this formula the gamma energy is < 1% of the total 
excitation energy of these events. The gamma energy correction is not included 
in calorimetry of this data as it is small compared to that of the excitation 
energy and is poorly constrained. 

3. Temperature determination with the momentum fluctuation ther- 
mometer 

A new method of calculating the system temperature can be derived from the 
momentum fluctuations of particles in the center of mass frame of the fragment- 
ing source. In this case, the frame is that of the reconstructed quasi-projectile. 
The momentum is constructed for each particle in the projectile frame center 
of mass for each axis. The momenta are then compared using Eq. |4] 



where P x , P y , and P z are the momenta in the X,Y and Z-axis respectively of 
each particle. The Q distributions extracted from equation 2] for the longitu- 
dinal (Q z ) and transverse (Q X y) are shown in figure O In an ideal, spherical, 
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Figure 2: Q z and Qx,y, Q distributions extracted from equation [4] for the longitudinal (Q z ) 
and transverse (Qxy) calculated using all fragments, (color online). 



fragmenting source the sum of Q for all fragments in the event in the center of 
mass of that event should be a S function at zero. However, in a real system 
fluctuations occur within a class of events. These fluctuations convert the mo- 
mentum S function into a distribution that can be characterized by a mean and 
width. In an equilibrated system, the mean should still equal zero. 

The variance (a 2 ) may be obtained from the Q z distribution through 

o 2 = (Ql) (Q z ) 2 (5) 

where Q 2 is the quadrupole moment (Eq.U]). If the mean equals zero, the second 
term cancels. Taking the first term as 

(Q 2 Z )= Jd* P (2P 2 -P 2 ) 2 f( P ) (6) 

and assuming a Maxwellian distribution of the momentum yields 

(Ql) = , ' 3/2 / d "P ^ - ^zPt + Pt) e-^^ 1 . (7) 
A Gaussian integral is used 

poo 

I n (a)= / x n e~ ax2 dx (8) 
Jo 

and the results that arc of interest to this derivation are: 

*>(«) = ^ W 
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It is important to note that these integrals are only for to oo and therefore 
should be doubled for -oo to oo. After derivation, the relation between variance 
and temperature becomes: 



where 



1 



(13) 



2mT 

For a single fragment type from a nuclear multi-fragmentation: 

a 2 = l2A 2 mlT 2 (14) 

where mo is the mass of a nucleon and A is the mass number of the fragment. 
For multiple particle types, the formula becomes: 

a 2 = \2mlT 2 Y,^A l ) 2 (15) 

i 

where £ is the concentration of the particle in question. A similar derivation 
can be carried out using Q x and Q y defined in Eq. @] 

Because of the velocity source cuts, the Z-axis does not exhibit the same 
momentum distributions as are seen in the X,Y-axes. Therefore, the temper- 
ature was re-derived using only the and P y momentum and employing the 
same Maxwellian distribution and Gaussian integral assumptions. Here Q 2 y is 
defined as 

Ql y = J d 2 p (Pi - 2P 2 P 2 + P$) f(p) (16) 

and using the Maxwellian distribution assumption gives: 

a 2 = Am 2 T 2 . (17) 

For a single fragment type from a nuclear multi-fragmentation: 

a 2 = AA 2 mlT 2 (18) 

where mo is the mass of a nucleon and A is the mass number of the fragment. 
For multiple fragment types, 

a 2 = AmlT^iCiAif (19) 

i 

where £ is again the concentration of a given particle type. The temperature is 
now linked to the variance through Eqs. [Inland |19l The temperature evaluated 
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Fi gure 3: Qx, Qyi Qz find Q^y (Ecjs. l4l and I18|l temperatures (Eqs. 1151 and I19|l calculated 
using all fragments as a function of E*/A of the source event, (color online). 



from all fragments assigned to the quasi-projectile fragmentation event is shown 
in Fig. |3] for the quadrupole oriented along various axes (Q z ,Qx,Qy) as well as 
the transverse quadrupole (Q X y) (Eq. |4j. These quantities should be exactly 
the same in a thermally equilibrated system, any deviation from thermal equi- 
librium will be disclosed from the differences in the quantity Q calculated along 
different axes. The statistical error bars are smaller than the size of the points 
in Fig. G3 The Q x , Q y thermometers are identical within errors and the Q z 
thermometer is slightly lower due to the longitudinal velocity cuts imposed as 
part of the quasiprojectile source selection. The transverse thermometer based 
on Qzj, provides a higher source temperature across all excitation energies than 
was obtained using the Q x , Q y , or Q z thermometers. This difference in tem- 
perature is the result of eliminating the effect due to the velocity cuts along the 
Z-axis. Hence it is imperative to investigate the temperature extracted from the 
transverse velocity of the fragments. 

The momentum distributions may also be widened by recoil during the 
source break up. The recoil effect on the momentum Q distributions has been 
parameterized as 

Qxy(0) — Qxy(app) I ^ _ ^ 

(20) 

where Q xy (o) is the recoil corrected value, Q xy (app) is the experimental value 
(Eq. [T5| , A is the mass of the source, and A/ is the mass of the fragment being 
considered fl3j |. 

Combining all of the fragments together to create a Q distribution depends 
on the assumption that all of the fragments, regardless of mass or charge, have 
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Figure 4: Apparent temperatures extracted from the Qxy distributions as a function of particle 
type, (color online) 



the same distribution width. This implies a simultaneous, statistical emission 
of fragments. The apparent temperatures obtained from the recoil corrected 
Q xy distributions as a function of particle type and source excitation energy 
are plotted in Fig. 01 Though the temperature increases with excitation energy 
regardless of particle type, there is a significant spread in measured temperature 
across the fragment types. 

A significant mass dependence at constant Z may be seen in the 3 ' He iso- 
topes. The 3 He yields a ~35% higher temperature than the 4 He across all E*/A. 
The 4 He apparent temperatures show a ~15% increase across E*/A= 1.5-7.5. 
The 3 He temperatures, on the other hand, increase by ~50% across the same 
E*/A region. It is possible that some of the mass dependence of the apparent 
temperature is the result of differential fragment emission time. It is known that 
3 He is emitted early in the fragmentation process. Conversely, 4 He is emitted 
throughout the de-excitation cascade p9j . The difference in measured temper- 
ature seen in the A=3 isotopes may indicate a Coulomb contribution to the 
momentum distributions widths. On the other hand the results for 2 H and 4 He 
isotopes rules out a strong Coulomb dependence of T. Other authors results [3] 
suggest that Coulomb is not so important in our beam energy region. 

In addition to the above mentioned factors, the temperature obtained with 
this thermometer could be a combination of thermal energy and Fermi momen- 
tum in the detected fragments [13j. Following the paper of Bauer, the measured 
fragments exhibit a momentum distribution resulting not only from thermal 
sources, but also from the Fermi momentum of the component nucleons with 
the fragment. This effect increases as a function of fragment size until a limiting 
value is observed. This correction is only meaningful for fragments with A>2. 
Thus, the momentum distributions of protons are free of this complication. 

The caloric curve derived from the proton momentum fluctuation widths of 
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the 86 Kr+ 64 Ni system is plotted in Fig. [3] (top panel). The fluctuations for 
this thermometer include the Z-axis, eq.(16), thus include collective effects if 
any. For reference, the compilation of Natowitz et al. [3] and two Fermi Gas 
curves ( T = v / E*/a with a = A/8, A/13 MeV -1 ) are plotted. As shown in 
this figure, the protons produce temperatures similar to what is obtained from 
a Fermi Gas at low E*/A indicating that this thermometer does indeed yield 
reasonable temperature values. The errors plotted in Fig. [5] for the Y-axis are 
estimated systematic errors of 0.5 MeV. 

As discussed above the quadrupole fluctuations displayed a nonequilibrium 
effect still present along the beam direction. We expect these effects to be 
present in the determination of the excitation energy as well. To a first ap- 
proximation we can correct for this nonequilibrium effect by replacing the ki- 
netic energy of the fragments in eq.(2) by its transverse component, i.e. : 
K cp — > 3/2Ktcp, similarly for neutrons. The excitation energy determined from 
the transverse kinetic energy of the fragments is approximately 0.3to 0.5 MeV 
larger than that which is determined by the total kinetic energy, see Fig. [5] 
(bottom panel). As we see in the figure, the difference when the transverse 
energy and the 'transverse' temperature are used becomes quite important for 
increasing excitation energies where collective effects become more important. 
Interestingly the 'caloric curve' obtained without collective effects shows a back 
bending similar to what was predicted by D. Gross long agof3lj. This might be 
a signature for a phase transition as discussed in [32l . |31| , however a more de- 
tailed study of this feature is beyond the goal of this paper and will be discussed 
in more detail in a future publication. 

When the collective motion is removed the transverse temperature rises from 
approximately 4 MeV to 10 MeV over this excitation energy range. Slope "tem- 
peratures" were also extracted from this data set by fitting the kinetic energy 
spectra[33] of the fragments. Both the quadrupole temperature and the slope 
temperature are derived from the kinetic energy. However the quadrupole tem- 
perature should be less susceptible to secondary decay and non-equilibrium 
effects. Since the velocity cuts imposed to minimize the effect of the non- 
equilibrated motion in the beam direction mainly affect the Z=l and Z=2 parti- 
cles, the slope temperatures extracted from the fits to the kinetic energy spectra 
for Z>3 should not be impacted. The values of the slope temperature were con- 
sistently approximately 3 MeV higher than those extracted from the quadrupole 
temperature. As already mentioned the fragment kinetic energy spectra are a 
result of both thermal energy and Fermi momentum [l3j]. Thus this compari- 
son may contain important information about the contribution from the Fermi 
momentum and will be the subject of future works. 

We expect that at low T or excitation energy the Fermi gas formula E* = 
aT 2 should be valid, see Fig. 5. Thus it is convenient to plot the quantity E* /aT 2 
versus T to see any deviation from a pure Fermi gas. In Fig (6) we plot the 
reduced excitation energy versus temperature as derived from the transverse 
kinetic energy of the fragments (full squares) as well as the total kinetic energy 
of the fragments (open squares). As seen in the figure at small T the two 
results are identical within the error bars. For increasing T the two results differ 
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Figure 5: Q z (Eq. |4j temperatures derived from proton momentum fluctuations as a func- 
tion of E*/A of the source event. For reference the caloric curve for A=60-100 from the 
Natowitz compilation [i| as well as two Fermi Gas ( T = yj aE* /A with a = 8, 13 ) curves 
are plotted(Top panel). (Bottom panel) same as above but for the corresponding 'transverse' 
quantities, (color online) 



much demonstrating the increased influence of the non-equilibrium effects with 
increasing T or excitation energy. The back bending is now more pronounced 
compared to fig. 5. 

The values obtained from the transverse energy have been normalized to 
1 at the lowest measured temperature using a level density parameter ao = 
A/13.3 MeV" 1 . We notice that such a quantity is density dependent i.e. a = 
ao(po/ p) 2//3 • From this relation, assuming that the nuclear Fermi gas is reaching 



different densities at different temperatures [34|, |35|, |36j , we can deduce a value 
of density for each measured temperature. We can easily calculate the density 
at each temperature by noticing that: 

— = (^l!)3/2 (21) 

Po V E* ' V ' 

The validity of the Fermi gas model is related to the ratio of the temperature 
to the Fermi energy T/e{p)p . If the density decreases also the Fermi energy 
decreases as: e{p)p — 36(p/ po) 2 ^ 3 (MeV). Where 36MeV is the Fermi energy 
of the nucleus ground state and we are assuming to be the same value at the 
smallest measured temperature. Since we can estimate the change of density 
for each temperature we can also obtain a value for the Fermi energy. In figure 
(7) we plot the reduced density versus the reduced T ( i.e. divided by the 
corresponding Fermi energy) . As we see in the figure the values of the reduced 
temperature are well below 1 which implies that our Fermi gas approximation is 
valid and the system is a quantum Fermi system even at the highest measured 
temperatures. This is at variance with many authors which assume that the 
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Figure 6: E*/aT 2 vs T derived from the transverse kinetic energy of the fragments (full 
squares) and the total kinetic energy (open squares) for the s6 Kr + 64 Ni system, (color 
online) 



system is classical at the highest measured temperature|30|. More important, 
the result without collective effects suggests that the density is roughly constant 
in the range of temperatures considered at variance with the analysis which 
includes collective effects. The latter case would suggest an expanding system 
at all temperatures while the previous case enforces the concept of freeze out 
volume. If this conjecture will be confirmed by future analysis we could study 
the nuclear equation of state at (almost) constant density. 

If the picture of an expanding system with increasing temperature is correct, 
a caloric curve as the one of figures (5) and (6) might be misleading. A different 
physical interpretation of the caloric curve could result if we plot opportunely 
normalized quantities as in figure(8). In these units the 'apparent' saturation 
of the temperature at high excitation energies now disappears because of the 
changing density of the system contained in the Fermi energy and the level 
density parameter. In contrast the results with non-equilibrium effects removed 
display a more remarked back bending. We stress that these results are derived 
assuming a Fermi gas relation and are for sake of discussion. A more precise 
determination of the equation of state requires a more complete determination 
of the quantities at play, i.e. excitation energy, temperature and density (or 
pressure) of the system. 

4. Concluding remarks 

Four momentum thermometers have been defined here (Eqs. [4j 1 19[) . The 
first three (Q x , Q,y, and Q z ) should theoretically be degenerate. However, the 
Z-axis was found to be truncated by the source selection. This affected all 
three of these thermometers. To remove this characteristic along the Z-axis, the 
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Figure 7: Reduced density vs reduced temperature derived for a Fermi gas from the transverse 
kinetic energy of the fragments (full squares) and the total kinetic energy (open squares) for 
the s6 Kr + 64 Ni system, (color online) 



temperature was then derived using only the X and Y-axes. This thermome- 
ter provides a slightly higher temperature reflecting the removal of the Z-axis 
with its incorrect width. The high value of T obtained from this thermometer 
may be a result of Fermi momentum in the fragments. This effect is zero for 
protons. The proton caloric curve is shown in Fig. [5] Since a non-equilibrium 
effect is present in the temperature determination, we have also suggested a 
method to reduce the non-equilibrium effects in the determination of the ex- 
citation energies by using the transverse kinetic energy of the fragments. The 
modified caloric curve thus obtained has been discussed in Fig. [5] demonstrating 
the increased importance of collective effects at high system temperatures (or 
excitation energies) . Assuming a Fermi gas relation for an expanding system we 
have been able to derive a density value for each temperature. From the density 
we can estimate the corresponding Fermi energy and demonstrated that the sys- 
tem is still a quantum Fermi system at the current temperatures and densities. 
We stress that the assumption of an expanding Fermi gas is suggested by some 
microscopic calculation |30j, but it is at variance with some statistical models 
which assume a constant freeze-out density [32| . If a constant density is the 
correct picture then the best representation of the caloric curve would be the 
ones of Fig. 5 and Fig. 6 while if density is changing we should refer to Fig. 8. 
In both cases we do not expect the densities to be dramatically different from 
those estimated in Fig. 7 which implies that we are in a quantum regime at all 
temperatures. 

Which picture of fragmentation is correct will probably be debated both 
experimentally and theoretically in future years. To this debate we have to 
add the observation that our proposed method gives different temperatures for 
different fragments. The origin of this difference could be in the statistics of 
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Figure 8: Caloric curve in reduced units for the °°Kr +° Ni system. Full squares are from 
the transverse kinetic energy of the fragments and the open squares are using the total kinetic 
energy, (color online) 

the fragments (i.e. Bose versus Fermi systems), in different Coulomb charges 
or time of creation including evaporation. Other proposed thermometers also 
suffer from this isotope dependence of the temperature. These features should 
be resolved in order to gain full understanding of an excited quantum nuclear 
system. 
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